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ABSTRACT 

We investigate the evolution in the distribution of surface brightness, as a function of size, for 
elliptical and SO galaxies in the two clusters RDCS J1252.9-2927, z = 1.237 and RX J0152.7-1357, 
z = 0.837. We use multi-color imaging with the Advanced Camera for Surveys on the Hubble Space 
Telescope to determine these sizes and surface brightnesses. Using three different estimates of the 
surface brightnesses, we find that we reliably estimate the surface brightness for the galaxies in our 
sample with a scatter of < 0.2 mag and with systematic shifts of < 0.05 mag. We construct samples 
of galaxies with early-type morphologies in both clusters. For each cluster, we use a magnitude limit 
in a band which closely corresponds to the rest-frame R, to magnitude limit of Mb = —18.8 at 
z = 0, and select only those galaxies within the color-magnitude sequence of the cluster or by using 
our spectroscopic redshifts. We measure evolution in the rest-frame B surface brightness, and find 
— 1.41 ±0.14 mag from the Coma cluster of galaxies for RDCS J1252.9-2927 and —0.90 ±0.12 mag of 
evolution for RX J0152.7-1357, or an average evolution of (—1.13± 0.15)z mag. Our statistical errors 
are dominated by the observed scatter in the size-surface brightness relation, a = 0.42 ± 0.05 mag for 
RX J0152.7-1357 and a = 0.76 ±0.10 mag for RDCS J1252.9-2927. We find no statistically significant 
evolution in this scatter, though an increase in the scatter could be expected. Overall, the pace of 
luminosity evolution we measure agrees with that of the Fundamental Plane of early-type galaxies, 
implying that the majority of massive early-type galaxies observed at z ~ 1 formed at high redshifts. 

Subject headings: galaxies: clusters: general — galaxies: elliptical and lenticular, cD — galaxies: 

evolution — galaxies: fundamental parameters — galaxies: photometry — galaxies: 
clusters: individual (RDCS J1252.9-2927, RX J1052.7-1357) 
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1. INTRODUCTION 

The apparent evolution in the colors and magnitudes 
of early-type cluster galaxies has long been used as 
a test of the history of galaxy evolution. The re¬ 
lations between magnitude and size or velocity dis¬ 
persion play an important role, as they give the ob¬ 
server a way to predict the apparent luminosity of 
a galaxy based on some other observable property 

jFaber/^^jarfeoi^^ Fisherirr977l: iKormendvl 

Il977t IDiorgovski^^Da^ 1198711 . Using these scaling 

relations, it has been observed that there is remark¬ 
able uniformity in cluster early-type galaxy properties 
at low redshifts, z 0-2. llSandag e fc Per elmuteilll991t 

l.lorgensen et al.ll996t iBernardi et al.l2003l for e^mmolel. 

At redshifts up to ~ 1, there have been a number of 
papers that have established a decrease in the mass-to- 
light ratio, using the Fundamental Plane. These studies 
find that the trend with redshift corresponds to a pas¬ 
sively evolving stellar population which has a luminosity- 
weighted las t epoch of major star-formation around 
z ~ 2 — 3 (Ivan Dokkurn fc FranxI 119961: iKelson et alJ 
I J 119981: IKelson J I2nnnii 

van Doki^im'^rTYanxnZflflli Iva.n DokTimT^Tstanfordi 

85721. 

W. M. Keck Observatory, 65-1120 Mamalahoa Hwy, Ka- 
muela, HI 96743 
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l2nn^ iWnvts et a,].l 12001 iHolden et a,l] l2?i?il . These 
systems also appear to have old stellar populations 
when examined using their colors ( e. q.\ |Bower_£LaD 
19921: lAragon-Salamanca. et alJ ITo^ lEllis et alJ llOOTt 
Stanford et al.ll998HBl^eslee et alJ2003bHHolden et all 

2004|i . with a similar last epoch of star formation of 
z ~ 2 — 3. Combining the evolution in colors and in 
the scaling relations provide us with the data to measure 
the luminosity-weighted age of early-type galaxies and, 
therefore, the epoch when most of the constituent stars 
formed. 

The size-magnitude relation is a projection of the 
Fundamental Plane (FP) along the velocity dispersion 
axis. It has long been known that l ow redshift early - 
type galaxies follow such a relation llKormendvl 1197711 . 
By examining how the apparent magnitude changes as 
a function of redshift at a fixed size, we can measure 
the apparent luminosity evolution in early-type galax- 
ies. Earlier measurements of this a t z ~ 0.4 — 0.6 
llSandag e fc Perelmuteill990albL 11 99111. and at z 

1.0 llSch ade et alJ 119961 I 199'A iSandage fc LubinI 

l200ltlLubin fc SandagelbooT m , found the same pace 
of evolution as found by the evolution of colors or from 
the FP of early-type galaxies. 

We investigate the size-surface brightness relation for 
two complementary goals. First, we can extend the 
measured evolution to z = 1.237, which corresponds to 
the h ighest redshifts that the FP has been measured t o 
date livan Dokkum fc Stanfordl2003HHolden et all200f)ll . 
However, because we have a much larger sample, we can 
also examine the distribution of the size-surface bright¬ 
ness relation to see if we find evidence for early-type 
galaxies that have younger stellar populations. Tradi¬ 
tionally, thi s has been done w ith colors or spectroscop¬ 
ically (see iTran et 120031 for “E-l-A” galaxies that 
have e arly-type morphologies). Ivan Dokkiim fc Stanford 
lf 2 f)nl . however, find that one of their early-type galaxies 
has an “E-l-A” like spectrum and appears to be an out¬ 
lier fro m the distribution of sizes and magnitudes for that 
cluster lIHolden et al.ll2?ioll . A significant population of 
such objects would imply that we are finding the epoch 
when intermediate mass cluster galaxies are completing 
their evolution into cluster early-type galaxies. 

Our sam ple has two cluster s. RDCS J1252.9-2927 at 
z = 1.237 (|RosatLet_aI] 1 20 041) an d RX J0152.7-1357 at 
z = 0.837 llDella Ceca et al.l200r)D . For each we use Hub¬ 
ble Space Telescope imaging data and other supporting 
ground-based data, discussed in § 2 , to construct a sam¬ 
ple of early-type cluster galaxies. For each galaxy in our 
sample, we measure three different total magnitudes and 
two different half-light radii (r^). The comparison of 
these measurements are discussed in §3. We verify that 
our sizes and total magnitudes are consistent, then use 
those two quantities to measure the mean surface bright¬ 
ness interior to the galaxy’s half-light radius (< fih >)■ 
In §4, we use r/, and < fih > io estimate the amount of lu¬ 
minosity evolution seen from low redshift to the clusters 
in our sample. In addition, we examine the scatter in the 
relation and search for galaxies out of the norm. We dis¬ 
cuss and summarize, in §5, our results in the context of 
other work. Throughout the paper we compute distances 
and absolute magnitudes using the cosmological parame¬ 
ters of ilm = 0.3, Ha = 0.7 and Hn = 70 km s~^ Mpc“^. 
We use the AH magnitude system lOkJlOOOD unless oth- 



Wavelength (A) 


Fig. 1.— A 13 Gyrs old early-type galaxy spectral energy distri¬ 
bution at a redshift oi z = 0.837, top, and 2 = 1.237, bottom. In 
both panels, we plot the filters we used for the observations with 
dotted lines and the Johnson B filter as a solid line in the rest- 
frame of the galaxy. In each case, we have a color that straddles 
the 4000 A break and a filter that overlaps with the Johnson B. 


erwise noted. 

2. DATA 

We combine a number of different datasets to ex¬ 
plore the evolution of the size-magnitude relation. For 
our high redshift samples, we use the Advanced Cam¬ 
era for Surveys (ACS) ima ging data for RDCS J1 252.9- 
2927, already discussed in iBlakeslee et al.1 ll2003bD . and 
for RX J0152 7-1357 . Both clusters are described in 
iPostma.n et al.l (j2nn5(l which focuses specifically on mor¬ 
phological classification and the morphology-density re¬ 
lation. Both clusters have additional ground-based data 
in the form of spectra and multi-wavelength imaging. We 
compare our high redshift re sult with the data fo r MS 
1358-1-62, at z = 0.328, from IKelson et all ll2000aD . We 
select this particular dataset because of the similarities 
in methods used by the authors. 

2 . 1 . RDCS J1252.9-2927 

As described in iBlakeslee et all ll2003tJl . the Advanced 
Camera for Surveys (ACS) imaged RDCS J1252.9-2927, 
at z = 1.237, for 12 orbits in the F775W filter, or 1775 , 
and 20 orbits in the F850LP filter, or Zgsp. The cluster 
was imaged as a mosaic of two by two positions, with 
three orbits in 1775 and five orbits in zgso at each posi¬ 
tion. The central area was covered by each of the four 
positions, resulting in deeper imaging over the central 
arcminute. In Figure ^ we plot the filter transmission 
curves and a representative spectral energy distribution. 
We note here that the 850 in the F850LP designation 
refers to the wavelength where the filter begins. The 
central wavelength of the filter, when combined with the 
CCD performance is 9077A, sampling at 4058A in the 
rest-frame of RDCS J1252.9-2927. 

The data were processed using APSIS, the ACS 
pipeline science investigation software, described in 
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iBlakeslee et alJ (l2003a^ . This pipeline generates drizzled 
images in each band, a detection image made by com¬ 
bining all of the drizzled data, and photometric catalogs 
constructed using a combination of the detection image 
and the final imag es in each filter. 

As discussed in IBlakeslee et al.l l|20n3bjl . we selected 
objects in RDCS J1252.9-2927 with zs 5 o < 24.8 to ensure 
adequate signal-to-noise for the fits. Eac h galaxy was 
fit with a single component Sersic model 
using vl.7a. Rev 3 of the GALFIT packa ge dPeng etlil] 
l2finl . The approach of iPeng et al.1 ll2002ll is to construct 
a model, convolve it with a point spread function and 
then compare the of the resulting model with the data. 
For a noise model, we use the variance maps produced 
by the APSIS pipeline. Our point spread functions are 
based on ACS imaging data of crowded star fields. We 
will discuss this further in ED 

Each galaxy with 0.5 < 1775 — zsso < 1-2 was mor¬ 
phologically identified by four of us (BH, NC, MP, and 
ME). These identifications w ill be discussed in another 
paper llPostman et al.ll200,5^ . The identifications were 
mad e using the T-type morphol ogical classification sys¬ 
tem llde Va.iicouleurs et al.ll199iTl . For the classifications, 
we found that 75% of the time all four classihers agreed 
on the morphological type, and three out of four agreed 
80% of the time. These statistics are for the whole of 
the sample and the agreement is much higher for the 
brighter galaxies. However, what matters for this work 
is how robust the classihcation for early-type galaxies in 
the final sample. There are only two galaxies in RDCS 
J1252.9-2927 that appear in our final sample where a 
majority of the four classifiers did not agree on a clas¬ 
sification as an early type (be.; lenticular or elliptical). 
Both of those galaxies appear near the magnitude limit 
of our sample and neither are part of the spectroscopic 
sample we will discuss below. At zsso < 23, 15 out of the 
20 early-type galaxies were classified unanimously. From 
24.5 > Z850 > 23, the fraction drops to 20 unanimous 
identifications out of 32 early-type galaxies, or 63%, but 
for 30 out of 32 times at least three out of the four clas¬ 
sifiers agreed, higher than the 80% average for the whole 
sample. 

Despite the high degree of internal consistency, one po¬ 
tential concern is a systematic offset in our classifications 
as compared with other publications. To control for this, 
one of the classifiers, MP, assign ed types to every galax y 
in the MS 1358-1-62 sample of iFabricant et ahl 11200011 . 
The classifications by M P of MS 1358-1-62 agree d 80% of 
the time with those from IFabricant et al.l 1200011 . 

In Figure [3 we plot the color-magnitude diagram of 
the early-type galaxies with T < 0 that are within 1 
Mpc, or 1192, of the cluster center and have zgso < 24.5 
mag. The circles represent galaxies classified as ellipticals 
(—5 <T< —4) and squares represent SO (—3 < T < 1) 
galaxies. The solid line is the color-magnitude relation 
{i—z) = 0.958—0.025(z85o— 23). The dotted diagonal line 
represents two s tandard deviati o ns, or 0.08 mag, from 
this relation, see IBlakeslee et all ll2003bll for details. We 
will use this color criteria for selecting early-type galaxies 
in the cluster. We remove spectroscopic non-members 
that, nonetheless, have the colors of the early-type red 
sequence in RDCS J1252.9-2927. We plot the magnitude 
limit for morphological identification, 2:850 = 24.5, as a 
vertical, dotted line in Figure |3 This magnitude limit 





Fig. 2.— Color-magnitude diagram for the ear ly-type galaxies 
in the ACS imaging data for RDCS J1252.9-2927 iBlakeslee et alJ 
l2QQ31Ji . The colors are measured within the half-light radius as 
measured from fitting Sersic models. The red circles represent el¬ 
liptical galaxies, while the blue squares are the galaxies classified as 
SO’s. The solid line is the mean color-magnitude relation while the 
parallel dotted lines represent the 2cr limits for the color-selected 
sample. The vertical dotted line is the magnitude limit for the 
color-selected sample. 


corresponds t o rs; 0.2£*. We calcul ate L* starting with 
the value from lNorberg et ^ 1200211 and evolvi ng it with 
the relation of lTOiT^oyamrfc Stanfo^ 1200311 . At 2 = 
1.24, M% = —21.7 in AB mag nitudes using the previous 
results ijBlakeslee et al.ir2003bll . 

In addition to the ACS imaging data described 
above, we also u se the photome tr ic and spectroscopi c 
data collected inIDemar col 1200,'^: IBosat i et al.l 1200411 : 
iLidman et al.l 1200411 : iToft et all 1200411 . The catalog 
for RDCS J1252.9-2927 contains photometry in the B, 
V, R, i, z, J and Kg bands. The near-i nfrared images 
were obtained with ISAAC on the VLT 
l2fiol . while the opt ical images were observed using the 
FORSl and FORS2 l|ADDenzeller fc RuDDrechtlll99^ in¬ 
strument on the VLT. These optical and near-infrared 
images were used to create a sample of galaxies for a 
redshift survey of potential cluster members. All galax- 
ies have Kg < 21 and J — K < 2.1 and R — K > 3 
ijDemarcd 120031 iRosati et ^1200511 . These colors were 
chosen to preferentially select galaxies at the redshift of 
the cluster regardless of the spectral energy distribution 
while ignoring redder stars as well as low redshift sys¬ 
tems. A total of 383 objects were observed, of which 235 
have secure redshifts. This yields a sample of 36 cluster 
members. The sample is almost complete for 2350 < 23, 
which is ~ 0.8L*. When comparing other samples with 
the redshift sample of RDCS J1252.9-2927, we will move 
the magnitude limit to 2850 < 23 but we will not use 
the color selection, shown with the two dotted lines, of 
Figure 13 This only adds one object to the sample that 
would not be included in the color-selected sample. 
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2.2. RX J0152.7-1357 

For RX J0152.7-1357, at z = 0.837, the ACS images 
were taken with the F625W, or rQ 25 , filter along with the 
1775 and zs 5 o filters. Each of the three filters was observed 
for two orbits. As with RDCS J1252.9-2927, RX J0152.7- 
1357 was observed using a two-by-two mosaic pattern 
with deeper imaging in the cluster core. In Figure 13 we 
show the inner 2' around the cluster center. This color 
image was made from the combination of all three filters. 

Our sample for RX J0152.7-1357 was constructed in 
as similar a manner as possible as our sample for RDCS 
J1252.9-2927. We used the SExtractor detection catalog 
to find galaxies. All objects with isophotal *775 < 24.8 
were fi t with a Sersic model with 1 < n < 4 and visually 
typed l|Postman et al.ll^005ll . The model fits were deter¬ 
mined in the *775 filter as that filter corresponds closest to 
the rest-frame B (see Eigure^, and the resulting aper¬ 
tures were used to measure the apparent colors. These 
results are shown in Figure 01 

For RX J0152.7-1357, we have obtained B, V, R, and 
I imaging data u sing the Keck t elescope with the LRIS 
ijOke et a'l.l1995|l . We used SOFI ) 

on the New Technology Telescope iDemarcoetall^OO^ ) 
to observe in the J a nd Ks bands. Using FO RSl and 
FORS2 on the VLT, iDemarco et al.l ll2rio5all targeted 
galaxies with i? < 24 using the photometric redshift 
selection criterion of 0.7 < Zp^ot < 0.95. 262 objects 
were observed resulting in 102 cluster members out of 
227 measured reds hifts. _ 

As discussed in iPostman et all ll2005ll . the classifica¬ 
tions for RX J0152.7-1357 were done by three of us (MP, 
NC, and BH). The error rates are effectively the same 
as for RDCS J1252.9-2927 but we examined all galaxies 
down to 1775 < 24, a fainter limit than our completeness 
for the redshift survey (discussed below). 

In the field of RX ,10152.7-1357. IDemarco et alJl|2005a|l 
found a group at z = 0.64 which has a significant number 
of early-type galaxies. These objects are, at best, slightly 
bluer than the cluster members (see Figure 0}). The in¬ 
clusion of these early-type galaxies in the photometric 
redshift selected sample reaffirms this. As such, we will 
not make a color selected sample for RX J0152.7-1357 
but, rather, will only use the spectroscopic redshift sam¬ 
ple. In the magnitude range of 20 < 1775 < 24, the com¬ 
pleteness of early-type galaxies ranges from 80% to 0%. 
The effective magnitude limit of the survey is 1775 < 23 as 
there is only one early-type galaxy fainter than that with 
a measured redshift. In the range 20.25 < ^775 < 22.75, 
roughly 70% of the early-type galaxies with spectra are 
cluster members. In that same magnitude range, 60% 
of all early-type galaxies have spectra. T his complete¬ 
ness is roughly the same as MS I358-I-62 l|Kelson et al.l 
l2flflflrjl against which we will compare our data. How¬ 
ever, it should be realized that the photometric redshift 
selection imposes, in effect, a color selection as photomet¬ 
ric redshifts are based on colors. Objects with stronger 
4000 A breaks are more likely to be selected for spec¬ 
troscopic follow-up. For the remainder of the paper, we 
select all early-type galaxies in RX J0152.7-1357 with 
1775 < 22.75 that are spectroscopically confirmed cluster 
members. This corresponds to ~ 0.5L*. 

2.3. MS 1358+62 


For our sample in RDCS J1252.9-2927 and RX 
J0152.7-1357, we used filters that straddled the 4000 A 
break for cluster members (see Figure^. To measure 
the amount of apparent luminosity evolution, at a fixed 
size, requires that we measure the magnitudes and sizes 
in filters that match, as closely as possible, the same 
rest-frame band pass. 

For a low redshift sample, we use the results fro m 
MS1358-I-62, at z = 0.33, from iKelson et al.1 l|2000albl ini 
because of the thorough study of the surface brightness 
profiles in the rest-fra me B band. For MS 1358-1-62, 
iKelson et al.1 (|2000albll conducted a spectroscopic survey 
to determine the velocity dispersion for a flux-limited 
sa mple of galax i es. T he 52 objects were selected from 
the lFisher et alJ lll998f) study of MS 1358+6 2 to R < 21. 
The total sample of IKelson et al.l l|2000albll includes an 
additional three galaxies that are below the magnitude 
limit, but we used only those with i? < 21. The selection 
for this sample was done regardless of morphology, so we 
include only galaxies that are classified earlier than 0, or 
the SO/a morphological clas s, according to the scheme of 
Ide Vaucouleurs et al idnii. 

We use the tabulated values of t he observed re s t-frame 
V band surface brightnesses from IKelson et aP ll2000all 
along with the half-light radii in that paper. To convert 
the V band surface brightness in to a B band surface 
brightness, we use the color from IKelson et all l|2000cll 
(taken from Ivan Dokkum et H] 1199811 . The apertures 
used to define these colors are the half-light radii mea¬ 
sured by fitting a de Vaucouleurs profile. In addition, 
the authors fit Sersic models. There is a difference in 
their Sersic models, however, w hen compared to the 
Sersic models in this paper. In IKelson et al.1 l|2000all 
the authors allowed the Sersic parameter to the range 
1 < n < 6. We restricted n over a smaller range, 
1 < n < 4 for our fits of RDCS J1252.9-2927 and RX 
J0152.7-1357. So, in all cases where lKelson et aP (|2000a|l 
finds n > 4, we replace the results of that fit, namely Vh 
and < uLh >, with corresponding n = 4 results from 
IKelson et all l|2000ajl. 

3. MEASURES OF THE SIZE AND MAGNITUDE 

In or der to measure th e evolution in the Kormendy re¬ 
lation ijKormeudvIll 97*^ , we need robust measures of the 
half-light radius and surface brightness. Below we discuss 
two methods, one parametric and one non-parametric, 
we used to measure both the total magnitudes, from 
which we infer the surface brightness, and the half-light 
radii. 

3.1. Model Fitting 

To determine the colors of each of the galaxies, we used 
models to establish the appropriate size of the apertures. 
We fit a single Sersic profile and allowed the parameter 
n to vary over the range 1 < n < 4. We imposed this 
restriction because we found that, if we did not, a few 
objects would settle to very high n numbers and, cor¬ 
respondingly large half-light radii. As the goal for the 
fits was to provide radii for measuring the colors of the 
red-sequence galaxies, large radii potentially could lead 
to systematic errors in the measurements. 

Th e models were fitted using GALFIT (|Peng et alJ 
|2002), and yield the half-light radius and total magni¬ 
tude of the model. The mean interior surface bright- 
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Fig. 3.— Color image made from combining the re 25 , *775 and 2350 imaging data for RX J0152.7-1357. The image is 2'on a side, with 
North up and East to the left. Some “red” galaxies in this image are actually members of a foreground group at 2 = 0.64. 


ness within the half-light radius for that model can then 
be computed as < Uh >= m + 5 log^g rh + 2.5 logj^g 27r 
ll.Torgensen et al.iri995ll . The GALFIT output for the size 
is actually the semi-major axis of the best htting ellip¬ 
tical model. We convert this radius into the effective 
radius, where a and b are the semi-major and 

semi-minor axes respectively, when we discuss the half- 
light radius. As mentioned above, the htting procedure 
constructs a model, convolves it with the point spread 
function, and then compares the resulting between 
the model and the data. Because the modeling process 
explicitly includes the point spread function, even for ob¬ 
jects with sizes smaller than the point spread function, a 
full width at half maximum of ~ 0V09, could be, in prin¬ 
ciple, reliably measured. However, errors in our point 
spread functions could cause systematic biases in our re¬ 
sulting half-light radii. Thus, we will exclude from our 
analysis objects with model half-light radii less than O'.'l 


for the rest of the paper. 

One of the advantages of the GALFIT software is the 
ability to ht multiple component models to multiple 
galaxies simultaneously. For most of the galaxies, how¬ 
ever, we masked out neighboring objects using the results 
of the SExtractor segmentation map, i. e., the map of 
pixels above the isophotal detection threshold. At min¬ 
imum, each postage stamp was 50 by 50 pixels. The 
actual size was four times the square root of SExtrac- 
tor’s isophotal area, the number of pixels one a above 
the sky. For close galaxies, however, we had to fit models 
simultaneously. This was often the case for the luminous 
ellipticals in the cores of the clusters. In cases where a 
multiple galaxies had to be fit, we ensured that all of 
the pixels that would be in the individual stamps using 
the above recipe were included in the stamp containing 
multiple galaxies. 

For every object in the color-selected and spectroscopic 
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Fig. 4 . — Color-magnitude diagram for the early-type galaxies in 
the ACS imaging data for RX J0152.7-1357. The colors are mea¬ 
sured within the half-light radius. Red circles represent elliptical 
galaxies while blue squares are the galaxies classified as SO. The 
dotted line represents the magnitude limit. Symbols with solid 
dots are confirmed to be cluster members while points with crosses 
are non-members. Unlike RDCS J1252.9-2927, we did not select 
galaxies using colors because of the presence of a 2 = 0.64 group 
of galaxies. Known members of this group can be seen as crosses 
with colors very close to the “red-sequence” of the RX J0152.7-1357 
members. 


redshift sample of RDCS J1252.9-2927 and spectroscopic 
redshift sample of RX J0152.7-1357 we fit a pure de Vau- 
couleurs, or n = 4 Sersic profile, in addition to the al¬ 
ready existing Sersic model. We selected a fixed n model, 
in addition to the Sersic fits, for two reasons. First, the 
models have one less degree of freedom and therefore 
the potential to be more robust. Second, a de Vau- 
couleurs’ model has traditionally been used. Thus, we 
can compare with a broader range of results from the 
literature. We plot in Figures [S] and 0 the observed size- 
surface brightness relation using the results from pure de 
Vaucouleurs model fits. The relation appears linear but 
with a large scatter for RDCS J1252.9-2927. This scatter 
likely comes from a combination of observational errors, 
the inclusion of non-cluster members and systematic bi¬ 
ases in addition to the intrinsic scatter in the relation. 
The size-surface brightness relation fo r RX J0152.7-1357 
is much tighter. The biweight scale l|Reers et a.1.ll199n|) 
of the scatter around the mean relation for RX J0152.7- 
1357 is 0.42 ± 0.05 for the de Vaucouleurs models and 
0.36±0.03 for the Sersic models as opposed to 0.76±0.10 
for the de Vaucouleurs and 0.72±0.09 for the Sersic mod¬ 
els for RDCS J1252.9-2927. We estir nate the errors oii 
the scatter using a iackknife process llBeers et al.lll990 l: 

iD^rr^ . 

In general, the error bars in Figures Eland El are some¬ 
what misleading. GALFIT reports the errors on the in¬ 
dividual terms but it is well known that there are sig¬ 
nificant anti-correlations b etween the errors for the sur¬ 
face brightness and size llHamabe fc Kormendvl ITo^ 
I.Torgensen et al.lli993H . Fortunately, these errors are al- 
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Fig. 5.— Size-surface brightness relation for the early-type galax¬ 
ies in the ACS imaging data of RDCS J1252.9-2927. Filled symbols 
are spectroscopically selected members and the open symbols are 
color selected, red circles for ellipticals and blue squares for SO’s. 
The surface brightnesses are in zs 50 magnitudes. We show the re¬ 
sults for Sersic models with 1 < n < 4 in the top diagram while 
the bottom diagram shows the sizes and magnitudes estimated us¬ 
ing de Vaucouleurs models. The best fitting relation to the whole 
sample is shown as a solid line, while the magnitude limit of the 
color-selected sample is shown as a dotted line and the magnitude 
limit we impose for the s pect roscopically selected sample is shown 
with a dashed line, see for details. We plot the semi-major 
axis of the ellipse of the anti-correlation between and < > as 

a line. The slope of the anti-correlation is very close to the slope 
of the size-surface brightness relation. We do note that the sample 
appears incomplete at the smallest sizes, <(7.'!. 


most parallel with the actual size-surface brightness rela¬ 
tion, so objects will be scattered along the relation. We 
illustrate the expected anti-correlation in Figures Eland 
El Both figures show a good deal of incompleteness for 
objects with sizes <071. For RX J0152.7-1357, this 
appears to be from the magnitude limit of the sample. 
However, for RDCS J1252.9-2927 we see a lack of objects 
with rh <07l and faint surface brightness even though 
such galaxies would be above our magnitude limit. We 
assume that this is a point spread function effect, as dis¬ 
cussed above, and will not use galaxies with sizes <07l 
for our analysis. 

3.2. Non-Parametric Magnitudes 

The main problem with both of our current measure¬ 
ments of the size and total magnitude is that they de¬ 
pend on fitting models to the data. To complement 
th ese model rnagnitu des, we imp l ement ed the approach 
of lWirth et al.l 11199411 and IWirt hI lll996ll. Thi^ approach 
uses the Petrosian rj measure llPetrosia3ll97fill to derive 
a model independent aperture. To review, the Petrosian 
rj = I{r)/ < I{< r) > computes the ratio of the surface 
brightness, I(r) at a particular radius, r to the average 
surface brightness within that radius < /(< r) >. This 
is the inverse of the function originally developed by Pet¬ 
rosian, but this version produces more tractable errors. 
Given rj for our galaxies, we need to decide on how to use 
























7 



j_jj _I_I_I_I_ .11.1 

0.1 1 


^half (arcsec) 

Fig. 6 . — Size-surface brightness relation for the early-type galax¬ 
ies in the ACS imaging data of RX J0152.7-1357, shown with filled 
circles. As with Figure[2| red circles are ellipticals and blue squares 
are SO’s. The surface brightnesses are in 2775 magnitudes. We 
plot the results for Sersic models with 1 < n < 4 in the top dia¬ 
gram while the bottom diagram shows the sizes and magnitudes 
estimated using de Vaucouleurs models. The solid line shows the 
mean relation while the dotted shows the magnitude limit. As 
with Figure we plot the semi-major axis of the ellipse of the 
anti-correlation between and < > as a line. RX J0152.7- 

1357appears less incomplete at small sizes, possibly because of the 
brighter magnitude limit imposed by the redshift survey used for 
selection. 


it to der ive an aperture. In the l iterature ther e are two 
choices. IWirth et alJ (119941) and IWirtH (119961) selected 
the aperture at which ?; = 0.1 to measure the magnitude. 
This aperture contains 95% of the total light of an expo- 
nen tial disk and 87% of the light of a de Vaucouleurs pro- 
file (IWir thlll996ll. The Sloan Digital Sky Survey (SDSS 
iStraiiss et al T ]2001 use Petrosian derived apertures as 
well, specifically the aperture containing twice the di¬ 
ameter at which r] = 0.2. This aperture contains 99% 
of the total light of an exponential disk but only 82% o f 
the light of a de Vaucouleurs profile (iStrauss et al l2n(Ul . 
Given that early-type galaxies are dominate d by de Vau¬ 
couleu rs profiles, we select the approach of IWirth et al.1 
119941) instead of th e approach of the SDSS outlined in 
IStrauss et al ] (l200l . 

When computing the n distribution , we do not fit el¬ 
liptical isophotes, as did IWirthI 119961) . but use circuly 
apertures as is done in t he SDSS IStrauss et ~ ] l200l . 
iLubin fc Sandagj l2001ar show that circular apertures, 
when plotted against effective radii, reproduce the re¬ 
sults of fitting elliptical isophotes. To compensate for the 
lower signal-to-noise at larger radii, instead of using the 
adaptive bin size approach of IWirthI 11996D . we choose 
to bin logarithmically in radius. Finally, to determine at 
what radius rj = 0.1, we interpolate our rj distributions 
for each object using a low order polynomial. 

This still leaves the question of correcting the aperture 
magnitude to a total magnitude. For each object, we 
measure the aperture where rj crosses 0.1 and then use 


that size to measure a total magnitude. This rj = 0.1 
magnitude is computed by interpolating from the curve 
of growth of the galaxy. This magnitude is not, however, 
the total magnitude of the galaxy. As we do not want to 
use a mo del-dependent aperture correction, we add 0.1 
mag to this magnitude. This would yield a magnitude 5% 
larger than the true total magnitude for an exponential 
disk and 3% less than the true total magnitude of a de 
Vaucouleurs profile. As most of our galaxies have Sersic 
values of 3 < n < 4, this aperture correction should, in 
principle, yield around 97% to 99% of the total light for 
the galaxies. For the rest of this paper, we will refer to 
this as our Petrosian total magnitude. 

3.3. Magnitude and Size Comparisons 

With two different models used for fitting and an addi¬ 
tional non-parametric total magnitude for each galaxy in 
each cluster, we can independently estimate the sizes of 
the errors on the sizes and total magnitudes. The statis¬ 
tical errors based on the expected variance in the photon 
flux for these galaxies are tiny, on the order of a hun¬ 
dredth of a mag. However, this ignores the errors caused 
by how well the model describes the data or the size of 
aperture. The formal errors in the fits are also small, 
with a P6r degree of freedom around 0.6 — 0.7 for the 
residuals. However, it is straightforward to construct 
models that fits a given distribution well, in the sense 
that the models yield a small per degree of freedom, 
but actually deviates systematically from the model. A 
good illustration of this are galaxies which are well fit by 
n = 4 models but, when fit with a Sersic, yield a value of 
n far from 4. Our non-parametric measure of the total 
magnitude checks for this. 

In Figured we show the distribution of the differences 
in magnitude for RX J0152.7-1357 as a function of the 
two different sizes. We show the difference between the 
total magnitude as measured by either the Sersic model 
(squares) or the Petrosian total magnitude (circles) and 
the de Vaucouleurs model plotted as a function of the 
half-light radius from the de Vaucouleurs model in the 
bottom panel. In the top panel, we plot the difference 
between the de Vaucouleurs model magnitude and the 
Sersic model magnitudes as squares while the solid dots 
are the difference between the Petrosian total magni¬ 
tudes and the Sersic model magnitudes but this time 
as a function of the Sersic model half light radius. We 
find 0.11 mag of scatter between the Sersic and the de 
Vaucouleurs magnitudes, 0.15 mag of scatter between the 
Petrosian and the de Vaucouleurs magnitudes and only 
0.06 mag of scatter in the differences between the Sersic 
and Petrosian magnitudes. In addition, there are system¬ 
atic shifts, the de Vaucouleurs magnitudes are on aver¬ 
age 0.02 ± 0.01 mag or 0.10 ± 0.02 mag brighter than the 
Sersic or Petrosian magnitudes respectively. However, 
the Petrosian magnitudes are only dimmer by 0.04±0.01 
mag than the Sersic magnitudes. We expect that the 
Petrosian total magnitudes should be somewhat smaller 
than the Sersic magnitudes by < 0.03 mag. 

We show the same quantities for RDCS J1252.9-2927 
in Figure IHl The scatter around the mean is larger, 0.17 
mag for the Sersic model magnitudes and 0.21 mag for 
the Petrosian magnitudes as compared with the de Vau¬ 
couleurs model magnitudes. As before, there are small 
systematic shifts. We find that the Sersic magnitudes are 
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Fig. 7.— The difference in the total magnitudes for RX J0152.7- 
1357 plotted as a function of size. The top panel plots the difference 
between the de Vaucouleurs model magnitudes and the Sersic mag¬ 
nitudes as squares, and the Petrosian minus Sersic differences as 
dots. This plot uses the Sersic half-light radius as the abscissa. 
In the bottom panel, we plot the difference between these Sersic 
model magnitudes and the de Vaucouleurs model magnitudes, as 
squares, and the difference between the Petrosian magnitudes and 
the de Vaucouleurs model magnitudes, as dots, as a function of the 
de Vaucouleurs model half-light radius. There ar e a f ew outliers 
when comparing with Petrosian magnitudes, see EH for details. 
Regardless of the magnitudes or sizes used, the overall scatter is 
small, ~ 0.1 mag, with no statistically significant offset between 
the different magnitude types. 
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Fig. 8 . — Same as Figure [3 but using the results of RDCS 
J1252.9-2927 at z = 1.237. Compared to RX J0152.7-1357, the 
scatter is larger, ~ 0.2 mag, and there appears to be a small cor¬ 
relation between the difference in magnitude and de Vaucouleurs 
model half-light radius. Most of this apparent correlation comes 
from galaxies with best fitting Sersic indices of n < 2 .5. This 
correlation, and the outliers, are more fully discussed in EH 



Fig. 9.— Four galaxies from RX J0152.7-1357. For each galaxy 
we show the original image (top), the residual image from the Sersic 
model fit (center) and the de Vaucouleurs model residuals (bot¬ 
tom). We give the Sersic index in the bottom left corner of the top 
image, and the A between the Sersic and de Vaucouleurs model 
magnitudes the top left corner. In some cases, there are still strong 
residuals after removing the Sersic model which either point to a 
higher Sersic index being a better fit, or galaxy substructure. 


0.03 ± 0.03 mag dimmer than the de Vaucouleurs mag¬ 
nitudes, while the Petrosian magnitudes are 0.07 ± 0.03 
mag dimmer. The scatter between the Sersic magnitude 
and the Petrosian magnitude is smaller, 0.16 mag, and 
the systematic shift is also smaller, 0.04 ± 0.02 mag. 

In Figures [71 and 0 there is an apparent correlation be¬ 
tween the plotted differences and some of the measured 
sizes. We examine the data using a Spearman Rank Cor¬ 
relation test. This test yields a parameter, referred here¬ 
after as t, which is the estimate of the degree of the cor¬ 
relation and can range — 1 < t < 1. Examining the data 
for RDCS J1252.9-2927, we find that t = 0.30 for the 
correlation between the de Vaucouleurs magnitudes mi¬ 
nus the Sersic model magnitudes as a function of the de 
Vaucouleurs half-light radii. The expected error on this 
measurement of the correlation, assuming that there is 
no true correlation for our sample is 0.09, yielding 3.3 
(T difference from 0 or no correlation. There is less of 
a correlation with Petrosian magnitude, t — 0.21 ± 0.09 
and we find no correlation, t = —0.05 ±0.09 between the 
difference in the Sersic and the Petrosian magnitudes as 
a function of the Sersic half-light radii. These positive 
correlations mean that the difference in the Sersic model 
and the de Vaucouleurs model magnitudes increases with 
large size, or the Sersic model magnitudes are fainter 
than the de Vaucouleurs model magnitudes as we mea¬ 
sure larger sizes. This correlation exists regardless of 
whether we use the Sersic model or de Vaucouleurs model 
half-light radii. These correlations do not appear for RX 
J0152.7-1357, however. For this cluster we find correla¬ 
tions with t = 0.12 ± 0.10 or smaller. 

Most of the apparent correlation in Figures 17181 comes 
from a handful of galaxies. In each case, the galaxy 
has a large de Vaucouleurs model magnitude when com¬ 
pared with the Sersic or Petrosian magnitude. We ex¬ 
amined each of these in detail. In every case where a 
galaxy lies more than three standard deviations from the 
mean of Figure jSl the best htting Sersic model has a 
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Fig. 10.— A histogram of the distribution of Sersic indices in 
RDCS J1252.9-2927, RX J0152.7-1357, and MS 1358-h62. There 
is, as expected, a strong peak at n = 4, as we constrain the fits 
to 1 < n < 4. There are a number of galaxies, however, that are 
classified as early-types but have small n values. Such galaxies are 
the strong outliers in Figure!^ 


low, n < 2.5 Sersic index and smaller half-light radius. 
Therefore, the de Vaucouleurs model is not a good fit 
to the data. In Figure [SI we show two examples of this, 
along with two examples of models where the Sersic and 
de Vaucouleurs models are in good agreement. In Fig¬ 
ure in we plot the distribution of Sersic indices for the 
three data sets, RDCS J1252.9-2927, RX J0152.7-1357, 
and MS 1358-1-62. In each case, the galaxies are identi¬ 
fied, by eye, as early-type galaxies. It appears that the 
higher redshift clusters have a broader range of Sersic in- 
dices than MS 1358-1- 62. This will be further discussed in 
iPostman et al.l ((200^. Nevertheless, the de Vaucouleurs 
models appear to be a bad fit to a number of galaxies, 
which should explain the scatter in Figure |H1 
When comparing the Petrosian magnitudes with the 
de Vaucouleurs model magnitudes for RDCS J1252.9- 
2927, there are four strong outliers. Each of these outliers 
is caused by neighboring galaxies that contaminate the 
Petrosian total magnitude. This is because our measure 
of the Petrosian total magnitude assumes nothing about 
the underlying surface brightness profile. We simply use 
the SExtractor segmentation maps to mask out bright 
pixels. We do not attempt to compute the contribution 
to the galaxy’s flux by the surrounding objects. Thus, 
in crowded areas where a larger number of pixels are 
contaminated by luminous or extended objects such as 
the brightest cluster members, our Petrosian estimates 
are significantly overestimated and so appear as large 
negative values in Figures 0 and |H1 
To compliment our Petrosian, or non-parametric total 
magnitudes, we can estimate a half-light radius by de¬ 
termining the radius where the curve of growth contains 
half of the total light from the galaxy. However, unless 
we make an attempt to remove the effect of the point 


spread function, we will And a bias towards larger sizes 
as the galaxies become intrinsically smaller. We com¬ 
puted these half-light radii and find that, on average for 
RDCS J1252.9-2927, the Petrosian sizes are 2% larger 
than the Sersic model sizes and 5% larger than the de 
Vaucouleurs model sizes for objects with model sizes of 
at least 0'.'15. This size bias grows markedly at smaller 
model radii, as expected. We find no statistically signif¬ 
icant difference between the de Vaucouleurs model sizes 
and the Sersic model sizes. For RX J0152.7-1357, we find 
very similar numbers. The scatter in the size as deter¬ 
mined by the de Vaucouleurs model as compared with 
the Sersic model is slightly larger because of the handful 
of cases where the best fitting Sersic parameter is low, 
n < 2.5 which causes large magnitude discrepancies as 
well. Overall, the agreement between the Petrosian es¬ 
timated half-light radii and the model radii are good. 
However, because of the bias towards larger sizes in our 
Petrosian half-light radii, we rely on the model half-light 
radii for the rest of this paper. 

Overall, it is encouraging that our different magnitude 
and size measurements yield a small scatter and little or 
no systematic offsets. Nonetheless, it is prudent not to 
mix the various measurements because of the small cor¬ 
relations seen above. The Petrosian magnitudes, which 
yield numbers close to the Sersic model magnitudes, 
could potentially be used in conjunction with the Sersic 
half-light radii. However, for the rest of the paper, we 
will measure evolution using the sizes and surface bright¬ 
nesses determined by the model fits and we will only 
compare Sersic model fits with other Sersic model fits, 
not with pure de Vaucouleurs models. 

3.4. Rest-frame Magnitudes 

To measure the amount of evolution, we need to com¬ 
pare the surface brightnesses in the same rest-frame 
pass band. Most measurements of the evolution of the 
FP are given in terms of the rest-frame Johnson B, or 
mass-to-light ratio. We used the templates of 
iBenftez et al ] il2nnl. which are modif ied versions of the 
templates from lCo1ema,n et af.l (|1 9S(Tl . to compute the 
Brest value as a function of an observed filter and a color. 
For RDCS J1252.9-2927, we computed Brest — -^sso as a 
function of *775 — ^850 while for RX J0152.7-1357 we com¬ 
puted Brest — *775 as a function of r625 — *775- For a given 
galaxy in RDCS J1252.9-2927 or RX J0152.7-1357, we 
then interpolate betw een the colors of the templates from 
IBenftez et all ll2004j) to then compute the Brest — zs 5 o or 
Brest — *775 color. Why we chose these particular filter 
combinations is illustrated in Figure 0 namely that the 
Z850 and 7775 filters match the Brest filter at the red- 
shifts of RDCS J1252.9-2927 and RX J0152.7-1357 re¬ 
spectively. The color was selected to span, roughly, the 
U — B color in the rest frame of the galaxy. This process 
is very similar to what was employed by iHolden et 
l|2nf)4ti . In that paper, the authors found that when the 
Brest filter was spanned by the obseryed colors, the er¬ 
ror on the magnitude transformation was on the order 
of 2%. In a case like RDCS JI252.9-2927, howeyer, the 
error was higher, at around 5%. 

4. APPARENT LUMINOSITY EVOLUTION 

4.1. Measuring the offset at a fixed size 
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We would like to investigate how much apparent evolu¬ 
tion there is in the size-surface brightness relation. First, 
we shall examine if the zero-point in this relation changes 
with redshift. The simplest way to measure this is to fit 
a linear relation to the size-surface brightness relation for 
the low redshift cluster, and then compute the constant 
required to minimize the deviation around the same slope 
for the higher redshift cluster. 

To determine the slope of the size-surface brightness 
relation at low redshift, we fit all of the data in MS 
1358-1-62. We then transform the MS 1358-1-62 data to 
the 1775 or zgso band in the manner described in E3I 
and rescale the sizes as appropriate for the z ~ 1 cluster, 
either RDCS J1252.9-2927 or RX J0152.7-1357. We limit 
both the z ~ 1 cluster and MS 1358-1-62 to the same mag¬ 
nitude. We then apply a series of cuts in the half-light 
radius, starting at a size of (y.'25 and working down to O'.'!, 
below which the effects of the point spread function will 
become important. At each size, we compute the zero- 
point in the size-surface brightness relation required to 
minimize the scatter around a line defined by the slope 
as measured for the whole of the MS 1358-1-62 sample. 
We compute this zero-point for both the trimmed MS 
1358-1-62 and z ~ 1 cluster samples. We then change the 
magnitude limit of the sample of MS 1358-1-62 by the 
currently measured amount of evolution. We trim the 
data at the new magnitude limit and continue iterating 
until we converge. To compute the zero-point for both 
distributions, we simply compute the biweight center for 
the deviations between the average surface brightnesses 
and the slope times the half-light radii. 

We find —1.04 ± 0.13 mag of evolution is required 
to match the size-surface brightness relation of MS 
1358-1-62 when comparing the de Vaucouleurs model re¬ 
sults for galaxies with rh >0'.'15 and zsso < 24.5 in 
RDCS J1252.9-2927. In rest-frame B, zgso = 24.5 cor¬ 
responds to Mb = —20.0 at z = 1.24 and Mb = 
— 18.8 at z = 0 after assumi ng the FP evolution of 
Ivan Dokkiim fc StanfordI ll2003r . We plot, in Figure ITTl 
the data from MS 1358-1-62 as open squares, both for 
the pure de Vaucouleurs models. The results for RDCS 
J1252.9-2927 are represented by solid dots. The solid 
line represents the fit to the whole data for MS 1358-1-62. 
The dotted line is the same but shifted by —1.04±0.13 B 
mag. This shift was calculated by minimizing the scatter 
around a line with the slope specified by the whole sam¬ 
ple. Repeating the process but using the Sersic model 
profiles, we find —1.08 ± 0.13 B mag of evolution (see 
Figure IT^. We note here, that, as we discussed above, 
iKelson et al.1 ll2000a|l allows a larger range of the param¬ 
eter n in their fits. All galax ies with n > 4 in Sersic 
models in iKelson et al.1 i2000all were replaced by the cor¬ 
responding de Vaucouleurs model. 

We make the same computation for RX J0152.7-1357 
and find —0.44 ± 0.06 mag of B band evolution for 
rh >0715, the size limit for the smallest galaxies in MS 
1358-1-62, and —0.66 ± 0.05 mag of evolution in the B 
band for the Sersic models. If we restrict ourselves to 
only those galaxies in RDCS J1252.9-2927 with redshift 
information and with zgso < 23, we find —1.41±0.16 mag 
(de Vaucouleurs model) and —1.42 ± 0.13 (Sersic model) 
mag of B band evolution for rh >07l5. 

We summarize the results of these fits in Tabled In 
that table, we also show the results for a selection of cut- 
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Fig. 11.— Size-surface brightness relation for the early-type 
galaxies in the ACS imaging data of both RDCS J1252.9-2927 and 
RX J0152.7-1357 using de Vaucouleurs fits. In each case, the low 
redshift data, MS 1358-1-62, are represented by the open squares 
while the filled circles are the data for the color selected sample 
from the redshift sample of RX J0152.7-1357 (top), the spectro¬ 
scopic redshift sample from RDCS J1252.9-2927 (middle) and the 
color selected sample from RDCS J1252.9-2927 (bottom). The 
galaxies in MS 1358-^62 are iteratively selected to be above the 
same magnitude limit plus the measured evolution to either RDCS 
J1252.9-2927 or RX J0152.7-1357. The solid lines represent the 
fits to the size-magnitude relation for all the data for MS 1358-1-62, 
while the dotted line is the relation shifted by the amount of mea¬ 
sured evolution. For RX J0152.7-1357, we find —0.44it 0.06 mag of 
evolution between 2 : = 0.328 and 2 : = 0.837. We find 1.41di0.16 mag 
of evolution for the redshift selected sample between MS 1358-t62 
and RDCS J1252.9-2927 while we find —1.04 di 0.13 mag of evolu¬ 
tion for the color-selected sample. 


off radii and for each of the samples of RDCS J1252.9- 
2927. We selected a number of half-light radii thresholds. 
We decided to not examine galaxies with rh <0710 as 
they are close in size to the full-width at half-maximum 
of the point spread function and the samples appear in¬ 
complete at those sizes, see EH However, there are 
few galaxies in the sample of MS 1358-1-62 that have 
rh <0715, so we selected that for our practical limit. 
We included a cutoff of rh >0720 along with rh >0710 
to illustrate that our results are insensitive to the exact 
threshold selected. 

The error we quote is the error on the mean shift be¬ 
tween the low redshift size-surface brightness relation 
and the high redshift relation. We assume no error on 
the slope measurement, simply taking the scatter around 
the mean relation and dividing it by the square root of 
the number of galaxies in our high redshift sample. This 
raises the question of how our particular choice of the 
slope of the size-surface brightness relation affects our re¬ 
sults. It appears, when looking at F iguresITTIandlT^ that 
the slope of the size-surface brightness relation is differ¬ 
ent for the different samples. We compare the data from 
MS 1358-b62 and RDCS J1252.9-2927 using the slope of 
the best fitting relation for RDCS J1252.9-2927 instead 
the slope of the relation for MS 1358-1-62. We find almost 
















11 



__l_I_ ■ .I_ l_ 

0.2 0.5 1 2 


(arcsec) 

Fig. 12.— Size-surface brightness relation for the early-type 
galaxies in the ACS imaging data of both RDCS J1252.9-2927 and 
RX J0152.7-1357 as compared with MS 1358-1-62. The sizes and 
sur face brightnesses are from the Sersic model profiles. See Figure 
m for an explanation of the symbols and the lines. The amount 
of measured evolution is similar as for the de Vaucouleurs model 
fits. We find —0.66 zb 0.05 mag for RX J0152.7-1357 (top) when 
the samples are compared with MS 1358-h62, —1.42 di 0.13 mag of 
evolution for the redshift-selected sample of RDCS J1252.9-2927 
(middle), and —1.08 di 0.13 mag of evolution for the color-selected 
sample of RDCS J1252.9-2927 (bottom). 


the same results with small changes, ~ 0.01 mag, in the 
measured evolution. For RX J0152.7-1357, our changes 
in measured evolution are larger, on the order of 0.03 
mag to, at most, 0.05 mag. This provides an additional 
source of uncertainty which we will add in quadrature. 

4.2. Comparison with the Fundamental Plane 

Ivan Dokkiim fc Stanford) ll2003ll measured the velocity 
dispersions for three galaxies in 

CL 08 48-1-4453, a low ma ss cluster of galaxies at z = 
1.276. IHolden et al.l II2005I1 measured velocity dispersions 
for four galaxi es in RDCS J1252.9-2927, whi ch, when 
combined with fvan Dokkum fc StanfordI ll2003ll . yields a 
sample of seven objects with z = 1.25. Using those val¬ 
ues along with sizes and surface brightness from HST 
imaging, the authors measured the offset in the FP from 
z = 0 to the redshift of the cluster. The evolution cor¬ 
responds to a change in the mass-to-light ratio (M/L) 
with redshift of —0.98z in the rest-frame B filter, or 
(—1.06±0.06) z mag of B band evolution at a fixed mass. 
For RDCS J1252.9-2927, this should correspond to —1.31 
mag of evolution in the rest-frame B band. To compare 
our relative measurements between clusters at different 
redshifts with those made using the FP, we need to com¬ 
pare our measured evolution with the Coma cluster , the 
z = 0 baseline for [van Dokkum fc Stanfor J ll2003ll and 
IHolden et al.1 ll200,5il" 

To measure how much evolution occurs between MS 
1358-1-62 and Coma, we implemented the same proce¬ 
dure as we used for comparing RDCS J1252.9-2927 and 


RX J0152.7-1357 with MS 1358-b62. For the data for 
the Coma cluster, we use the B half-lig ht radii and sur¬ 
face b right nesses, listed i n Table A 7 of lJorgensen et all 
lll996tl . In iKelson et alJ ll2000all . the authors compare 
the techniques used to measure the half-light radii and 
surface brightnesses in IJorgensen et al.1 lll996ll with the 
two-dimensional fitting techniques they used and find no 
significant offsets. Thus, we compare the Coma radii and 
surface brigh t nesses directly with the values tabulated in 
IKelson et alJ ll2000dl . We measure a shift of—0.35±0.07 
mag in the rest-frame B band between Coma and MS 
1358±62. These measurements were done with the whole 
sample of Coma and MS 1358±62, as both are magnitude 
limited, redshift selected samples. 

The amount of evolution expected between z = 0 and 
z = 0.328 is —0.35 ± 0.06 m ag in rest-frame B using 
the b est fitting relation from Ivan Dokkum fc StanfordI 
([2f)nl . in good agreement with the —0.31 ± 0.12 mag 
expected from Ivan Dokkum fc Franx and the 

—0.33 ± 0.12 mag from IKelson e^T i 199711 . Our result 
of —0.35 ± 0.07 mag is in excellent agreement but with 
smaller statistical error bars. The sy s temati c errors are 
the same as the FP, see IKelson et al.l Il2000all . 

We add the measured evolution between Coma 
and MS 1358±62 to the measured evolution we find 
for RDCS J1252.9-2927 and RX J0152.7-1357. We 
plot both in Figure 1131 along wit h th e FP results 
from Ivan Dokkum fc StanfordI ll2003ll and IHolden et alJ 
(1200,511 . We fit a line to the da.ta for RX J0152.7-1357 
and for the color-selected sample of RDCS J1252.9-2927 
and find (—1.13 ± 0.15)z mag of B evolution, very close 
to the (—1.06 ± 0.06)z mag found with the FP. We find 
(—1.42 ± 0.33)z mag of evolution if we use the redshift- 
selected sample of RDCS J1252.9-2927, a result not sta¬ 
tistically different from the FP because of the much error 
estimate. 

Overall, our agreement with the measurements from 
the FP are good, with the caveat that it is clear that 
sample selection plays a critical role. Our first sample 
of early-type galaxies in RDCS J1252.9-2927 was con¬ 
structed by finding all galaxies with the correct mor¬ 
phologies that also fall within the color-magnitude se¬ 
quence of RDCS J1252.9-2927. Restricting our sample 
to galaxies with sizes of Vh >0715, we find a result in al¬ 
most perfe ct agreement with restuls from the Fundamen- 
tal Plane (Ivan Dokkum fc StanfordI l200.3t IHolden et alJ 
[2f)f)l . regardless of whether we use Sersic or de Vau¬ 
couleurs models. If, however, we restrict ourselves to 
those galaxies identified only by redshift, we find more 
evolution at the highest redshifts. The redshift selec¬ 
tion was made using a K magnitude-limited sample, with 
some color selection. As shown in Figures|21and[Sl almost 
all early-type galaxies with z < 23 are included. There¬ 
fore, it is puzzling why such a sample would appear to 
show so much more evolution. 

4.3. Scatter in the Size-Surface Brightness Relation 

In Figure d we plot the the histogram of the scat¬ 
ter around the size-surface brightness relation. This is 
the difference, in magnitudes, between the measured sur¬ 
face brightnesses and those predicted by the size-surface 
brightness relation. We plot the residuals for RDCS 
J1252.9-2927 as compared with those from MS 1358±62, 
using the de Vaucouleurs model sizes and surface bright- 
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Fig. 13.— Amount of B luminosity evolution from the change 
in the zero-point of the size-surface brightness relation with red- 
shift. The results for MS 1358-1-62, represented as an open blue di¬ 
amond at 2 = 0.33, shows the measured amount of evolution with 
respect to Coma using our size-surface brightness measurements. 
The orange square represents the average of the Sersic and de Vau- 
couleurs results for the r^, <0^.^15 sample of RX J0152.7-1357. The 
solid red circle is average for the spectroscopically-selected sample 
of RDCS J1252.9-2927 and the open red circle represents average 
for the color-selected sample. We plot, with a solid line, a fit to 
the evolution of MS 1358-^62, RX J0152.7-1357, and the color- 
selected sample of RDCS J1252.9-2927. The dotted line shows 
the fit to the the same data except using the spectroscopically- 
selected sample of RDCS J1252.9-2927. For both lines, the er¬ 
ror bars at 2 = 1 represents the error in the slope. The results 
for either fit agree with the FP measurements in CL 0848-1-4453 
dvan Dokkiim fc StanfordI 120031) ■ show^ with a green cross, and 
RDCS ,71252.9-2927 liHolden^ al.l200Bli . shown with a purple tri¬ 
angle. 


nesses. The solid line represents the results for RDCS 
J1252.9-2927, while the dotted line is scatter for the com¬ 
parison sample of MS 1358-1-62. We plot the histogram 
containing only conhrmed cluster members in the top 
panel while the bottom panel shows all color selected 
galaxies. In both panels, the histograms are normalized 
to one. For this hgure, we plot all galaxies with pro¬ 
jected sizes of Th > 0'.'15 at z = 1.237. The magnitudes 
for the galaxies in MS 1358-1-62 have been brightened by 
the measured —1.04 ± 0.13 mag of evolution in the B 
band for the bottom panel and —1.44± 0.15 mag in the 
top panel. 

Examining Figure0again, however, we see that the ob¬ 
jects below the mean relation have, on average, a larger 
error in both the size and average surface brightness than 
those above the mean relation. Thus, a low tail towards 
large deviations is expected. At brighter magnitudes 
there are three galaxies brighter than expected based on 
the distribution of MS 1358-1-62. The observed scatter 
is, as mentioned before, 0.76 ± 0.10 mag around a fixed 
size. Given that the errors in our total magnitude mea¬ 
surements are, based on the comparisons between the de 
Vaucouleurs model magnitudes and the Petrosian total 
magnitudes, around 0.2 mag, most of the scatter must 
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A/j, (rest-frame B) 

Fig. 14.— The deviations from the mean size-surface brightness 
relation for RDCS J1252.9-2927 (solid line) and MS 1358-1-62 (dot¬ 
ted line) plotted as normalized histograms. All galaxies in RDCS 
J1252.9-2927 with >0^15 are included. The top panel contains 
all galaxies in the redshift sample and the bottom panel contains 
all those in the color selected sample, all galaxies within the dotted 
lines in Figure|^ The scatter for RDCS J1252.9-2927 is 0.76±0.10 
mag for the color selected sample and 0.65 ± 0.12 mag for the red- 
shift selected sample, in contrast to the 0.57 ± 0.07 mag found for 
MS 1358±62. 


be intrinsic or systematic. For RX J0152.7-1357, there 
are 55 early-type galaxies cluster members that meet our 
magnitude and size criteria. In Figure lTHl we plot the dis¬ 
tribution of the two sets of galaxies, and note that the 
scatter is 0.42±0.05 mag. For comparison, we find in the 
sample from MS 1358±62a scatter of 0.57 ± 0.07 mag. 

We calculate for each of our scatter measurements an 
er ror based on the ja ckknife of the data as recommended 
bv iBeers et al.l (1199011 . It seems that the scatter measure¬ 
ments are statistically very different. For example, the 
difference between the scatter in RX J0152.7-1357 and in 
MS 1358±62 is 0.15 ±0.06 mag, a 2.5 standard deviation 
difference. If we make the assumption that scatter in the 
size-surface brightness relation comes from a Gaussian 
parent population, we can use the F-ratio test to com¬ 
pare the variances of two different scatter measurements. 
When we compare the scatter for RDCS J1252.9-2927, 
0.75 ± 0.10 mag, with that of MS 1358±62, 0.57 ± 0.07 
mag, we find that, despite the 0.19 mag difference, the 
odds of the two populations being drawn from the same 
parent population is around 2%, not a statistically sig¬ 
nificant difference. We find a better agreement when 
comparing RX J0152.7-1357 with MS 1358±62, a 15% 
probability that variance measurements of the two sam¬ 
ples are drawn from the same distribution, despite the 
apparent 2.5 standard deviation difference. 

In H2.ll we discussed the accuracy of our morphologi¬ 
cal classifications. At fainter magnitudes, we found that 
our scatter in classification was higher, as measured by 
the disagreement amongst the classifiers. The inclusion 
of later-type galaxies could cause a bias in the estimate 
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Fig. 15.— The deviations around the mean size-surface bright¬ 
ness relation for RX J0152.7-1357 (solid line) and MS 1358-1-62 
(dotted line). This Figure is similar to Figure m All galax¬ 
ies with redshifts in RX J0152.7-1357, solid points in Figure 
and rh >0^.M5 are included. The scatter for RX J0152.7-1357, 
0.42 di 0.05 mag, is smaller than the results for MS 1358-^62. 


of the evolution. We measured the evolution using only 
those galaxies with unanimous agreement amongst the 
four classifiers and found a larger result, —1.46±0.13 mag 
for the color selected sample. This agrees well with the 
—1.44± 0.15 mag found for the redshift selected sample. 
The inclusion of late-type galaxies should bias towards 
more evolution, as spiral structure and large disks usu¬ 
ally hold star formation, yet the inclusion of the more 
uncertain classifications actually lessens the amount of 
evolution. However, the lower surface brightness galaxies 
are also the galaxies with more uncertainty in their clas¬ 
sification. Thus, by removing the galaxies with higher 
uncertainy in their classification, we naturally remove 
the lower surface brightness galaxies. The scatter for the 
color-selected sample, with uncertain objects removed, 
is 0.64 ± 0.12 mag, very close to 0.65 ± 0.12 mag we 
find for the spectroscopically selected sample. Thus, we 
raise the possiblity that the difference in the two sam¬ 
ples could come from classification errors as the majority 
of uncertain classifications are in the 23 < zgso < 24.5 
mag range, and below the cutoff for the spectroscopically 
selected sample. Removing these less certain classifica¬ 
tions lower the scatter and raises the amount of evolution 
measured in the color-selected sample. However, these 
objects have, on average, a lower surface brightness, 0.35 
mag in zssoi than the objects with unanimous classifica¬ 
tions. These lower surface brightnesses could also explain 
why the classifications are uncertain. 

5. DISCUSSION 

We observe, using the size-surface brightness relation, 
(—1.13±0.15)zi3 mag of luminosity evolution from z = 0 
to z = 1.237. If we restrict our highest redshift sam¬ 
ple to only those galaxies with spectroscopic redshifts, 
instead of the large sample of galaxies with red colors. 


we find (—1.42 ± O.SSjz mag of evolution, which is in 
statistical agreement with our other result. If we use 
only the color selection for RDCS J1252.9-2927, we find 
that both of the measurement s agree well with the FP 
results at the same redshifts (Sii.JDQkkum et al.lll998t 


Ivan Dokkum fc Stanfordl l2003l: lHolden et al.ll 2005() and, 

thus, with the overall trend in luminosity ev olution found 
from other size-surface brightness st udies llSchade et alJ 

ment points towards the conclusion that the majority of 
massive elliptical galaxies at these redshifts are passively 
evolving, in agreement with their color evolution. 

Given that we find slow luminosity evolution for 
the majority of cluster early-type galaxies in our sam¬ 
ples, this is another indication that the luminosity- 
weight ed_redslnf^_of_Jormntion is hig h for these ob¬ 
jects. Ivan Dokkum fc Stanford! ll2003ll finds Zf > 2.5, 
though for a differ ent cosmology, using the FP, while 
iHolden et al.1 l|200,5|l finds Zf2.3t.Q'^ for certain assump¬ 
tions about early-type galaxy evolution. As the evolu¬ 
tion found in that paper agrees with the evolution in 
offset of the size-surface brightness relation, we would 
find a similar r esult. This epoch of for mation is also in 
agreement with iBlakeslee et al.l (I2003IJ| from the colors 
of galaxies in RDCS J1252.9-2927. In general, the large 
scatter on the size-surface brightness relation prevents us 
from making significantly more accurate measurements 
oi Zf. Since we selected galaxies for red colors in RDCS 
J1252.9-2927, we could very well be missing an important 
part of the cluster early-type population, namely those 
galaxies that will become red early-types but are bluer 
than our selection criteria at the epoch of the observa¬ 
tions. This would bias our results to higher redshifts. 

The redshift selected sample of RDCS J1252.9-2927, 
however, shows much more evolution, (—1.42 ± O.SSjz 
mag as compared with the color-selected sample. It is 
possible that this is a selection effect. Namely, it is eas¬ 
ier to measure redshifts for higher surface brightnesses 
objects. We attempted to mitigate this by trimming the 
redshift-selected sample at zs 5 o < 23, where the redshift- 
selected sample is almost complete, containing 17 out of 
20 early-type galaxies with rh >0715, regardless of color. 
One possible explanation for this result is the scatter 
around the mean size-surface brightness relation. 

We observe that the apparent scatter for RDCS 
J1252.9-2927 around the size-surface brightness relation 
is larger than for RX J0152.7-1357. This result is not 
statistically siginificant though an increase in the scat¬ 
ter would be expected at higher redshifts. There are a 
couple of lines of evidence that, at z ~ 1, there are ellipti¬ 
cal galaxies that show evidence of recent star formation. 
Such objects would show a sm a ller th an av erage M/L. 
Both Ivan Dokkum fc StanfordI ll2003ll and iTran et all 
mm find galaxies morphologically classified as early- 
type galaxies but with strong Balmer absorption lines 
indicating recent starbursts. T he example galaxy in the 
Funda mental Plane sample of Ivan Dokkum fc StanfordI 
((2f)f)l has a much lower than expected M/L value, a 
factor of :;:^ 13 lower than the m assive galaxies in that 
sample, and IHolden et all (120051) finds a large scatter in 
the M/L value for a purely luminosity selected sample, 
much larger than the scatter for lower redshift clusters 
such as Coma or MS 1358-1-62. The inclusion of a handful 
of such galaxies in our sample would increase the scatter 
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but disentangling them requires spectra or good infor¬ 
mation. Spectroscopic evidence for recent bursts of star- 
formation in RDCS J1252.9-2927, which could explain 
both the larger scatter and the larger apparent evolution 
when using only the spectroscopic ally-selected samp le, 
will be discussed in a future paper l|Rosati et al.ll2001^ . 

6 . SUMMARY 

We have investigated the size-surface brightness re¬ 
lation for two clusters of galaxies, RDCS J1252.9-2927 
at z = 1.237 and RX J0152.7-1357 at z = 0.837. For 
each cluster, we identified, by eye, a sample of early-type 
galaxies across a broa d range in color and in magnitude 
llPostman et al ] |20^ . Using our ACS imaging, which 
covers the rest-frame B, we created a magnitude limited 
sample for each cluster. The magnitude limits corre¬ 
spond to approximately 0.3L* for RDCS J1252.9-2927 
and 0.5L* for RX J0152.7-1357. For RX J0152.7-1357, 
we sel ected galaxies b ased on the 102 spectroscopic red- 
shifts llDemaTcnll2flfl3^. RDCS J1252 . 9-2927 has 36 spec- 
trosc opic members (iDemarcd 120031: iRosati et all l2004 
l 2 ?iol and we expanded the sample for that cluster 
using a color selection. Our color criteria for RDCS 
J1252.9-292 7 utilizes the color-mag nitude relation and 
scatter from iBlakeslee et ^ (l2003blf . including all galax¬ 
ies within two standard deviations of the mean relation. 
In each case, we chose the observed magnitude limit in 
the filter that corresponds most closely to the rest-frame 
B band. 

We fit both Sersic and de Vaucouleurs models to each 
galaxy in our two cluster samples. These models pro¬ 
vide total magnitudes and half-light radii with the point 
spread function broadening removed. We converted the 
total magnitudes to the average surface brightness within 
the half-light radius. In addition to the model total mag¬ 
nitudes, we implemented a third total magnitude mea¬ 
surement. We used the P etrosian magnitudes as im¬ 
plemented in iWir IB iTT^ and found excellent agree¬ 
ment between these non-parametric measurements and 
the Sersic model magnitudes. The disagreement with 
the de Vaucouleurs models in general involved a small 
systematic of 0.04 — 0.1 mag depending on the cluster. 
In addition, we found that the difference between the 
de Vaucouleurs model magnitudes and either the Sersic 
or Petrosian magnitudes correlated with the model sizes 
for one of the two clusters. In general, we find that our 
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overall trends are robust, regardless of which size and 
magnitude we used for determining the surface bright¬ 
ness as long as we compare a given cluster measurement 
with a measurement made using the same model in a 
different cluster. 

We measured the amount of rest-frame B magnitude 
evolution by computing the shift in apparent surface 
brightness at a fixed size. We compare the c lusters in this 
paper t o the results from MS 1358-1-62 in IK el son (d 
ii2nnnaii . which contains both Sersic and de Vaucouleurs 
average surface brightnesses and half-light radii. For MS 
1358-1-62, we compute the amount of evolution required 
to match the size-surface brightness relation between MS 
1358-b62 and either RDCS J1252.9-2927 or RX J0152.7- 
1357, with the data trimmed at a fixed size and mag¬ 
nitude limit. We apply this correction to the data in 
MS 1358-1-62, recompute which galaxies would fall above 
our magnitude limit, and repeat until we converge. We 
also apply thi s to MS 1358-1-62 and the B band Coma 
cluster data of l.lorgensen et ^ l|199611 . For that, we mea¬ 
sure —0.35 ± 0.07 mag of evolution, in excellent agree- 
ment with the Fundamental Plane, or FP, measurement s 
l|van Dokkum fc Franxlll99'^ iKelson et al.lll99'^ l2000cjl . 

The rest of the evolution measurements, summarized in 
Table □ show results consistent with the FP in most 
cases. We find that the average surface brightness at a 
fixed sized changes with redshift as (1.13±0.15)z, very 
close to the (—1.15±0.13)z found with the FP. Only when 
we examine the subset of galaxies in RDCS J1252.9-2927 
selected by spectroscopic membership do we find a devi¬ 
ation. 

The overall trend in the luminosity evolution that we 
find agrees with this. At the redshift of RDCS J1252.9- 
2927, we are likely within a 2 —3 Gyrs of observing epoch 
when most of the stars formed in cluster early-type galax¬ 
ies. 
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TABLE 1 

Summary of Measured Evolution from MS 1358+62 


Cluster 

z 

Memberships 

Min. 

A/i (n = 4)'' 
{Brest mag) 

A/i(l < n < 4)'' 
(Brest mag) 

RX J0152.7-1357 

0.837 

z 

0'.'20 

-0.42 ± 

0.06 

-0.67 ± 0.06 

RX J0152.7-1357 

0.837 

z 

0'.'15 

-0.44 ± 

0.06 

-0.66 ± 0.06 

RX J0152.7-1357 

0.837 

z 

O'.'IO 

-0.57 ± 

0.06 

-0.68 ± 0.05 

RDCS J1252.9-2927 

1.237 

color 

0'.'20 

-0.97 ± 

0.14 

-1.17 ± 0.13 

RDCS J1252.9-2927 

1.237 

color 

0'.'15 

-1.04 ± 

0.13 

-1.08 ± 0.13 

RDCS J1252.9-2927 

1.237 

color 

O'.'IO 

-1.09 ± 

0.11 

-1.06 ± 0.11 

RDCS J1252.9-2927 

1.237 

z 

0'.'20 

-1.32 ± 

0.16 

-1.28 ± 0.11 

RDCS J1252.9-2927 

1.237 

z 

O'.'15 

-1.41 ± 

0.16 

-1.42 ± 0.13 

RDCS J1252.9-2927 

1.237 

z 

O'.'IO 

-1.53 ± 

0.16 

-1.47 ± 0.13 


^This column notes whether the galaxies in the sample were selected based on spectroscopic 
redshifts (z) or based on the galaxies’ colors. 

'^Minimum half-light radius for inclusion in the sample. 

^Average change in surface brightness measured for the sample based on de Vaucouleurs 
models (n = 4) or Sersic models (1 < n < 4.) 



